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We describe the properties of the neutron rich nuclei around N=28 in the shell mode frame- 
work. The valence space comprises the sd shell for protons an the pf shell for neutrons with- 
out any restriction. Good agreement is found with the available experimental data. The 
N=28 shell closure, even if eroded due to the large neutron excess, persists. The calcu- 
lations predict that '"^S and ''^S are deformed with f3 — 0.29 and P = 0.32 respectively. 

I. INTRODUCTION 

The exploration of the behaviour of the nucleus under extreme conditions; high spin, vicinity to the drip lines, finite 
temperature, etc. is a major source of new insights in nuclear structure. One of the most important questions raised 
by the study of the drip lines is whether the basic shell order evolves with the neutron (proton) excess An example of 
this problem is given by the N=20 very neutron rich nuclei where the strong anomalies detected were explained 
as a collapse of the standard N=20 shell closure. Extended shell model calculations |^-^ showed that the inversion of 
the spherical closed-shell configurations and the 2p-2h collective intruders could account for the large binding energies 
and half-lives, unexpected ground-state spins, and high level densities at low excitation energies. Also, mean-field 
calculations using the Skyrme interaction predicted a sharp transition from spherical to prolate shapes |^ 

Recently, there has been an increase of interest in the N=28 isotones far from the stability, motivated by; i) 
the possible existence of anomalies in the shell closures, as already found in N=20 and ii) the unusual '*^Ca/'*^Ca 
abundance ratio measured in the solar system or the Ca-Ti-Cr anomalies observed in some meteoritic inclusions. 

Sorlin et al. ^ undertook the study of the /3-decay and /3-delayed neutron emission probabilities for "''^S and 
45-47(^j^ Their measured half-lives were much larger than those predicted by TDA or QRPA calculations. 
These discrepancies were attributed to unexpected shape transitions in the region. Mass formulas such as the finite- 
range liquid drop , the finite-range droplet or the Extended Thomas-Fermi with shell corrections predict 
the existence of shape coexistence around ^^S. In all the cases the energy surfaces are very soft with close lying minima 
corresponding to different deformations. 

Similar conclusions have been drawn by Werner et al . Relativistic mean field calculations and non-relativistic 
ones using the Skyrme forces produce very flat energy surfaces for the Sulphur isotopes, with several minima separated 
by energy barriers of just a few hundred keV. In this his situation it seems necessary to go beyond the mean-field 
using, for instance, the generator coordinate method. 

In this paper we discuss all the neutron rich nuclei with N> 20 and 14 <Z< 20 in the Shell Model framework. 
This is, to some extent, complementary to the other methods mentioned above and can give new information on 
the structure of these nuclei. The paper is organized as follows. Section II gives a brief description of the model, 
discussing the choice of the valence space, the effective interaction and the coulomb corrections. Section III compares 
our results with the experimental data. In section IV we examine the behaviour of the N=28 shell closure far from 
the stability. Section V deals with deformation in heavy Sulphur isotopes. 

II. THE MODEL 
A. The Valence Space 

It consists of the full sd shell for Z-8 protons and the full pf shell for N-20 neutrons. Although the calculations reach 
very large dimensionalities, they can be performed without truncations in all cases. We do not include the intruder 
states that could be formed by exciting particles from the sd to the pf shell. These kind of states, that appear at 
low excitation energy around ^°Ca, should be less important near N=28. The reason is twofold; sd protons get more 
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and more bound with increasing neutron excess then unfavouring proton intruders; and the correlation energy of the 
configurations with neutron intruders do not gain much energy compared to the standard ones (the opposite is true 
in e.g. ^^Mg). This is in agreement with the observed spectra. The existence of many collective low lying states near 
N=Z=20 contrast with their absence for N~ 28, Z> 20. 



B. The interaction 



The construction of the effective hamiltonian is critical in any shell model description. It is well known that 
the saturation properties of the realistic interactions must be corrected phenomelogically in order to have a good 
description of the nuclear structure [ p^ . Unfortunately, the experimental data for the very neutron rich nuclei with 
14 <Z< 18 are scarce. Even if the masses are known over a relatively wide region, at N=28 they have only been 
measured for nuclei with Z> 19. As an example the masses of the Sulfur isotopes are only known up to N=26. Half- 
lives are only known for ^''S and ''^"'''^Cl. The level schemes of Ar and CI with 20 < < 23 are partly established. 
The situation is clearly better for the K isotopes where the evolution of of the 3/2+, 1/2+ doublet with the neutron 
number is known up to N=28 [ p^ . 

Our starting interaction has three parts. The USD interaction of Wildenthal [|l^ is used for the particles lying in 
the sd-shell. A modified version of the Kuo-Brown interaction, denoted KB in ref |po| , gives the two body matrix 
elements for the pf shell particles. This interaction is better suited to this valence space that the KB3 interaction 
that has been used in conventional p/ shell calculations |^^. Furthermore, when used in the lower part of the pf shell 
it also gives good spectroscopic results. Finally, the cross-shell interaction is the G-matrix of Lee, Kahanna and Scott 
(LKS) [Q. This is the part of the interaction whose monopoles we have proceeded to modify phenomenologically. 
The data used in the fit are: 



• The natural parity states of nuclei with few pf particles and sd holes. 

• The evolution of the 3/2+, 1/2+ doublet along the K isotopes. 

• The position of the 5/2+ states in "'^K. 

The experimental data cited above are not enough to fix all the monopole parameters. They essentially fix the 
interaction among neutrons in the lf7/2 and 2p3/2 orbits with the sd protons. On the contrary, the interaction with 
lf5/2 and 2pi/2 particles is not well determined. In our choice of the final interaction '^^Ca and '^^Si have very similar 
single neutron spectra as predicted in Dufio and Zuker ||2^ mass formula. There exists other possible choices of this 
part of the interaction that induce important modifications in the structure of heavy Silicons and to a minor extent 
in that of the heavy Sulphurs. 

The 2-body matrix elements incorporate a tiny mass dependence and are written as: 



Lo{Ao) 



< r2(A) > 



(1) 



where Aq = 40 and huJa=ll MeV are used to calculate the LKS G-matrix. Instead of the usual law < >= roA^^^ 
fm, we borrow from reference ||2^ a modified expression, which is better adapted to large mass regions. 



< >= To 



A 1-C 



N - Z 
A 



l/2^ 



-,1/3 



(2) 



leading to the following expression for the 2-body matrix elements 



V = 



1/2N 



-I 1/3 



Vn 



(3) 



with C = 0.42 and Aq = 40. 
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C. Coulomb Energies 



The absolute energies predicted with the precedent interaction are relative to an ^^O core and do not include the 
coulomb energy of the protons. Thus, to compare our predictions with the experiment, they must be corrected with 
the following prescription: 

E,iN, Z)=Ea + E,{N, Z) + Ef^'{N, Z) (4) 

where i^o is the inert core energy (127.62 MeV), Ec{N,Z) is the coulomb energy relative to the core and E^^ is 
the nuclear energy obtained with our interaction. The Ec energy, which is assumed to depend only on N and Z, has 
been calculated with the Antony and Pape formula P3| 



F(A 7)^ + 1.44(Z - l/2)Mi/3 _ 1.02 Z > Z, 

'^'^ ' ^ EciA,Z + l)-l.U{Z+l/2)/A^/^ + 1.02 Z < Z, ^ ' 



where Zg = A/2 or {A + l)/2 and the coulomb correction Ec{Zs) along the stability line is given by 

p (7 )- 0.162^2 + 0.95Z - 18.25 Z < 20 , . 

J^ci^s)- 0.125^2 + 2.35Z- 31.53 Z > 20 

which agrees with the Chung- Wildenthal coulomb correction ||^ when Z^ < 20 and with that of reference |Q for 
Z, > 20. 

III. RESULTS AND DISCUSSION 
A. Comparison with the data. 

In this section we compare compare the theoretical predictions with the experimental data. The shell model matrices 



are dealt with by the code ANTOINE |25| a very fast and efficient implementation of the Lanczos algorithm . The 
USD predictions for the sd-shell nuclei have been largely reported However, we shall include sd nuclei in the 

S2N plots to stress the continuity of the shell model description across the N=20 magic number. 

Table |, contains the predicted and experimental energies of the 2+ states, that show an excellent agreement. The 
interaction, adjusted to reproduce the behaviour of the 3/2+, 1/2+ levels alongthe odd-A K, predicts for the doublet 
a similar behaviour in the neutron rich Chlorine and Phosphorous (see table ||). The crossing of these two levels is 
due to the gradual degeneracy of the 2si/2 and ld3/2 orbitals as the occupation of the lf7/2 shell increases. In the 
P isotopes, the 1^5/2 orbit is nearly closed, consequently the evolution of the doublet gives us an idea of the relative 
position of the two orbitals. As we shall comment later, their relative possition is determinant for the development of 
deformation in the heavy Sulphurs. 

The spectroscopic information on the nuclei with Z< 20 N> 20 is limited to a some K, Ar and CI isotopes. Figures 
|l| to show the theoretical and experimental spectra for the cases -40-42j^^ 39-41^^. ^^^^ ^^C\- for which the energy 
levels are known. Yet, there are still many uncertainties in the measured level schemes, specially the spin assignments. 
The theory-experiment agreement is perfect for 40-4ij^ ""^Ar. The few states without theoretical counterpart are 
most probably sd-intruders. For the other nuclei the concordance is also quite good. 

In addition to these spectra, the half-lives of several Ar and S isotopes are known. Their decays are dominated by 
allowed Gamow Teller transitions. The ground states of the father nuclei were calculated in our Ohco space and the 
Gamow- Teller sum rule state was determined in a larger space that includes all the Ip — Ih in order to recover the 
3(A^ — Z) sum rule. The strength functions were obtained by the Whitehead method The phase space factors 



of Wilkinson ||27[] were used to compute the half-lifes. Total strengths and half- lives are compiled in table III. The 
agreement is very satisfactory. The calculated half-life of "''^S is very sensitive to the Qp. We have used the theoretical 
value; if we had employed the extrapolated number given in the Audi and Wapstra compilation ps| ], the half-life 
would have increased to 150ms. 

The two neutron separation energies S2Ar are known over a large region, unfortunately, for most Z's they do not 
reach N=28. The agreement with the experimental results is good as it can be seen in figures @ to |l2|. However, the 



discrepancies with the extrapolated masses |29 are important 
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B. Does the N=28 shell gap persist? 



To answer to this question we analyse our predictions for the S2Ar, the excitation energies of the first 2+ states, and 
the occupation numbers. Structural changes, such as spherical to deformed shape transitions, should show up in the 
behaviour of these quantities as a function of the neutron number. 

The S2Ar's around N=28 have some common features for different Z- values. In three cases: Ca, K and Si, a 
sharp decrease of S2Ar is observed at N=28. For the other isotopes, there are not sharp slope changes, but S2Ar still 
decreases. These results are consistent with a shell closure. This conclusion is opposite to that coming from the mass 
systematics, that predicts almost constant S2A' around N=28, and thar could be interpreted as a signature of a new 
region of deformation. 

Figure |l^ shows the evolution of the 2+ energy for the Ca and Si isotopes. The large energy increase of the 2+ 
for the four nuclei ^°Ca, '*^Ca, ^''S and ^^Si is something that typifies doubly magic nuclei. Notice, however, that the 
effect is smaller for the exotic ^^Si. A similar evidence comes from the lf7/2 occupation numbers along N=28. Those, 
shown in figure reveal an increase of the f7/2 occupation at Z=20 and Z=14, again slightly in "^^Si. 

The 2+ excitation energies of the isotopes of Ar and S, around N=28, are presented fig. [ij and compared with the 
corresponding energies of nuclei much closer to the stability at those neutron numbers (Ti, Cr). There is an upward 
shift of the excitation energies at N=28 that is a little smaller in Ar and S. The behaviour is seen to be very similar 
close and far off stability, thus indicating that the N=28 closure stands in the very neutron rich region. 

From the analysis of these three sets of results: S2Ar, E2;(2"'") and the occupation numbers we can conclude that the 
lf^^2 neutron configuration at N=28, although somewhat eroded by the large neutron excess, is still dominant and 
that ^^Si is a new doubly magic nuclei. Therefore, the vanishing of the shell closure far from stability that happens 
at N=20, does not seem to occur at N=28. 

C. The deformed nuclei and "^S 

Our calculations show the existence of low lying 2+ states for the Ar and S isotopes with 22 <N< 26. However, 
only the Sulphur isotopes have large (and negative) spectroscopic quadrupole moments and enhanced E2 transitions 
(see figs. 18 and 19). ^^Ar and ^'^Ar have reasonably large BE2's but very small spectroscopic quadrupole moments. 

The lowest 2"*" energies correspond to ^"S and ^^S. Their spectroscopic quadrupole moments are -IT.le^fm^ and 
-19.2e^fm^ respectively. Assuming that the Bohr-Mottelson formula for the intrinsic quadrupole moment is applicable, 
we get Qo=59.8 e^fm^ for '^"S and Qo=67.2 e^fm^ for ^^S, corresponding to deformation parameters /3 ~ 0.29 and (3 ~ 
0.32. Our results are consistent with the Relativistic Mean Field calculations of Werner et ai, which predict f3 ~ 0.25. 
Their non-relativistic calculations point to a more complicated situation with nearly degenerate prolate (/3 ~ 0.25) 
and oblate {f3 ~ —0.24) solutions. The only indication of coexistence in our results is provided by the presence of a 
very low (1.51 MeV) second 0+ state in '^^S. 

The values of the E2 reduced transition probabilities are 77.9 e^fm^ (A=40) and 93.0 e^fm** (A=42). We compare 
the quadrupole moments of ^^S and other nuclei near N=28 in figure |l4[ The value for "^s is about 2/3 of that 
predicted for ^*Cr leading to very similar deformation parameters (3 =0.32 and /? =0.30. In addition, the yrast levels 
of "^^S follow the rotational law J(J+1) up to J=6. 

It is also important to notice that, although a truncated lf7/2 2p3/2 valence space for the neutrons, contains the basic 
degrees of freedom for the description of deformed nuclei in this region, it enhances either the rotational behaviour 
or the closed shell features, depending on the basic configurations that are involved. As an example let consider the 
evolution of the 2+ excitation energies of the Sulphur isotopes. In the truncated space the energy of the 2+ in N=26 
goes down 200 keV to 0.88 MeV while in N=28 it goes up 300 keV to 1.78 MeV. 

In order to explain the onset of deformation we shall explore which are the configurations that are responsible for 
the quadrupole coherence in our space. It is well known that the SU(3) symmetry is broken by the strong I ■ s coupling 
in the pf shell. Yet, it has been recently shown that the subspaces of a major shell spanned by Aj = 2 orbits 
develop an approximate SU(3) symmetry, named quasi-SU(3). For instance the lf7/2 2p3/2 valence space is enough 
to to explain the collective behaviour of nuclei in the pf shell. In this space the maximal quadrupole collectivity is 
obtained for configurations with 4 protons and 4 neutrons. In our case, the basic neutron degrees of freedom are 
contained in this space, but protons are confined in the sd shell. For the nuclei we are interested in, the orbit ld5/2 
is essentially closed and only the ld3/2 and the 2si/2 can be considered active. In this situation neither SU(3) nor 
quasi-SU(3) can develop. 

However, one must look more carefully at the behaviour of these two orbits. As far as the neutrons are filling the 
sd shell, both orbits remain well separated. For instance, the excitation energy of the 3/2i level in ^^P and "^^P are 
1.39 and 2.62 MeV respectively. Thus, the ground state wave functions are dominated by configurations with the 
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maximum allowed number of particles in the 2si/2 orbit. Particularly, for ^^S the (2si/2)^ configuration represents 
the 89% of the wave function. When the neutrons begin to fill the pf shell, the two orbitals become more and more 
degenerate (see table |) inducing strong mixing in the wave functions. A valence space with degenerate ld3/2 and 
2si/2 orbitals has the geometry of pseudo-SU(3), and maximal quadrupole coherence is obtained for two protons, i.e, 
for the Sulphur isotopes. For three particles the quadrupole moment is strictly zero and beyond it changes of sign 
leading to oblate shapes. Thus, we expect other nuclei as, for instance, the Ar isotopes to have small colectivity. 

Notice, however, that strict degeneration is achieved only at N=28, where the neutron collectivity is rather small. 
In other words, there is a mismatch between the region where protons have collective wave functions (N=28) and the 
region were neutrons have quadrupole collectivity, that prevents the existence of a larger region of deformation. A 
compromise seems to take place for Z=16 and N=26 (''^S) or N=24 ('"'S). 



IV. SUMMARY 



In this paper we have studied the N=28 region far from the stability. We have built an effective interaction 
for sd active protons and pf active neutrons that is the natural generalization of successful sd-only and p/-only 
interactions. The agreement of the theoretical results and the experimental data is very satisfactory. However, much 
more experimental information is needed to obtain a complete understanding of the region. Accurate information on 
the 2+ excitation energies, on the evolution of the 3/2+, 1/2+ doublet along the CI and P chains, the single particle 
spectrum of "^^Si are challenges for the physics of exotic nuclei. We have also addressed the question of the magicity 
of N=28. A detailed analysis of the S2n, 2+ energies and occupation numbers allow us to conclude that it persists in 
the very neutron rich regime, contrary to what happens at N=20. Finally our results show that the ^o^zg isotopes 
are well deformed, but the effect seems to be very local; no large region of deformation is therefore expected. 
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TABLE I. 2^ excitation energies (in MeV). 



N 22 24 26 28 30 



Or 


tti u.oby v.liz u.yyj i.ozi u.oyy 
exp 0.869 0.751 0.783 1.434 0.835 


1 I 


m 1.209 0.869 1.003 1.755 1.103 
exp 1.083 0.889 0.983 1.555 1.047 


L^a 


-t-h 1 KQA 1 /inn 1 QQO Q O/IT^ 1 IIH 

til i.Or54 1 .4Uy i.ooo o.y4 < 1. 1 iU 

exp 1.524 1.157 1.347 3.832 1.030 


Ar 


th 1.368 1.249 1.182 1.655 1.219 

(.rp 1.4(il 1.208 - - - 


S 


th 1.209 1.047 1.073 1.465 1.272 
exp — — — — — 


Si 


th 1.900 1.670 1.652 2.558 1.461 
exp — — — — — 



TABLE 11. Energy difference between the 3/2f and l/2f states in K, CI and P. 



\N 20 22 24 



26 



28 



30 



K 

mm -mm 

Cl 

E{l/2) - E{3/2) 
P 

mm -mm 



th 2.73 0.96 0.67 

exp 2.52 0.98 0.56 

th 1.82 0.71 0.06 

exp 1.73 0.40 - 

th 2.61 1.32 0.68 

exp — — — 



0.28 
0.47 
-0.21 



-0.38 
-0.36 
-0.05 



0.19 



-0.06 



0.32 0.06 0.38 



TABLE III. Total Gamow- Teller strength and theoretical vs. measured half-lives of the even-even Ar and S isotopes. 



GT~ T*'' jiexp 

2 2 

44 Ar 22.85 4.70m 11.87m 

28.84 6.94s 8.4s 

16.78 174.4m 170.3m 

22.87 6.67s 8.8s 

28.89 389ms 560ms 

34.93 54ms 123ms 
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FIG. 1. Theoretical and experimental natural parity states for '"^A' 
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FIG. 3. Same as figure § for ^'^K 
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FIG. 5. Same as figure |l| for ^°Ar 
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FIG. 6. Same as figure |l| for *^Ar 
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FIG. 7. Same as figure § for ^^Cl 
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FIG. 9. S2n energies for the K isotopes. 
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FIG. 10. S2n energies for the Ar isotopes. The experimental points at the vertical line beyond are not measured quantities 
but numbers from Audi and Wapstra's extrapolations 
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FIG. 11. S2n energies for the CI isotopes. See caption to figure 
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FIG. 13. S2n energies for the P isotopes. See caption to figure 
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FIG. 15. 2+ excitation energies for the Ca and Si isotopes. 
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FIG. 17. £7/2 occupation numbers. 
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FIG. 18. Spectroscopic Quadrupole moments for the Cr,Ti,Ar and S isotopes. 
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FIG. 19. BE2 for Ca, Ar, S and Si isotopes. 
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